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c I '. Abstract 

C$ , We describe the asymptotic behavior of a filtration problem from a contaminated 

porous medium to a non-contaminated porous medium through thin vertical fissures 
of fixed height h > 0, of random thinness of order e and which are e-periodically 
distributed. We compute the limit velocity of the flow and the limit flux of pollutant 
at the interfaces between the two porous media and the intermediate one. 

o ■ 

OO 

1 Introduction 

We consider a porous medium which is contaminated by some pollutant and which com- 
municates with another porous and non-contaminated medium through vertical fissures 
of height h > 0, of random thinness of order e > 0, and which are periodically disposed. 
Each of these porous media is given a e-periodic structure. Instead of considering a Stokes 
problem in each of these porous media, we assume that the velocity of the fluid is governed 
by a Darcy law with periodic permeability matrix. 

The purpose of this work is to determine the influence of the fissures on the transport 
of the contaminant into the non- contaminated medium, computing the global flux of 
pollutant which penetrates in this non-contaminated medium and the asymptotic velocity 
of the fluid which flows through the fissures. 

Let O be a bounded, smooth and open subset of M 3 , with boundary T, such that 

o+ = nn {x 3 > o} 0, 
n h = Q n {x 3 < -h} ^ 0. 



*This work has been supported by the Comite Mixte Franco-Marocain under the PHC MA/08/183, 
t Corresponding author. Tel.: +33 3 89 33 63 10; fax: +33 3 89 33 63 19. 
Email addresses: Alain.Brillard@uha.fr (A. Brillard), m.eljarroudi@uae.ma (M. El Jarroudi). 
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Let E x {0} = dQ n {X3 = 0}. E is a bounded and smooth subset of R 2 . We define 



r r 





r+ 
r- 



Ex{0}, 
E x {-h}, 

{x e£l\ (xi, x 2 ) £ S, — /i < a; 3 < 0} , 

dQ+\r+, 

The domain Q is thus equal to Q + U Tq" U Y^ U IT U f^. 

Let (II, T, P) be some probability space and (T(t)) tgR be a group of transformations 
on (II, T), that is satisfying 

T(0) = Id n , 
T(tx+t 2 ) = T(tx)oT(t 2 ) Vh,t 2 eR, 
P {T- 1 (t) A) = P (A) WA £ X, Vi £ R, 

where Idn is the identity map on II and the set {(t, u) £ R x II | T (t) u G A} is dt x dP 
measurable, for every A £ T. We suppose that T is ergodic (or metrically transitive), 
which means that every A £ T such that T (t) A = A, for every t £ R, has a probability 
P(A) equal to or 1. 

We introduce some random processes q and r defined on IR x IT and satisfying the 
following conditions: 

1. q(t,u) is a stationary random process, that is, for every positive integer n, every 
points £].,••• , i n , every i in R, and every B e B (R), one has 

P ({w | g (t + t 1; w ) , • • • , q (t + t n , u) £ B}) 

= P{{" I q(ti,T(t)co),--- ,q(t n ,T(t)u) £ B}) , 



where i3 (R) is the Borel cx-algebra on R. Since T preserves the measure P, the 
above equality implies that the joint distribution of {q (ti) , ■ • • ,q (t n )} is the same 
as the joint distribution of {q (ti + t) , • • • , q (t n + t)} , for every t in R. 

2. The derivatives ^£ and ^£ exist for m = 1,2,3 and there exist non-random con- 
stants ci, C2 and C3 such that the following bounds hold true with probability 1 



< ci < g (t, u) < c 2 < 1 ; |r (£, w) | < 1 



d m q 



dt r ' 



d m r 



dt r 



< c 3 , 



(1) 



From the properties of T and q, we derive the ergodic theorem (see, for example, [TU] 
and [18]) 



\/n £ IT : (g" (0)) = lim — / g" (t, w) (it, 
T^ooZI J _j- 



T 



(2) 



almost surely, where the symbol (.) stands for the mathematical expectation with respect 
to the measure P. 
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Let (ctj (w)) ieZ and (/^ (w)) i€Z be sequences of random variables satisfying 

|a»H| <c 4 ; <c 4 , Vi eZ, (3) 

with probability 1, where c 4 is a non-random constant. We define, for every i,j e Z, the 
fissure F £i jj (w) as 



ea 



i (-£ e x 3 ) < x 1 



is < ea] 



(-e- e x 3 ) 



(-e e x 3 ) < x 2 - j'e < eat (_ e ^ 3 ) , x 3 e 0[ J ' 



with < £ < 1, < < 2/3 and of (2) = r (z + & (w) , w) ± g (z + a t (u) , u) /2. Let 
J £ (w) = G Z 2 I F £j jj (a;) C Y^}. We also define the sets 



0r e>ij (w) n r + , 



a £ H = <9F e (w) \ (r+ £ u r~ £ ) (w) n(w) 



u rJ £; ..H, 

(z,j)e/ e (w) 
(jj)e/ £ (aj) 

(i,j)6J e (w) 



Let Z = ]-l/2, l/2[' 3 be the unit cube of M 3 and assume that it can be decomposed as 
Z = Z 1 U S U Z 2 , where Z 1 and Z 2 are two disjoint, open and connected sets separated 
by the smooth surface S (Fig. 1). 



Figure 1: A 2D view of the periodic structure of the porous media. 

We assume that Q + and are two e-periodic porous media which communicate 
through the fissures Y £jij (cu). 



Figure 2: The porous media and the fissures. 



3 



We set 





= n+n 


u 


< w 


= Q+n 


U 
V fcez 3 




= Q + n 


U 





= flfc 


n 


u 

V fcez 3 




= n- h 


n 


U 


s; 


= flfc 


n 


u 

V fcez 3 



We suppose that Vt^ ,e (resp. fi^'J) is the portion of f2 + (resp. Q h ) consisting of the 
pores which are filled in with some fluid and Qp £ (resp. fi^f) is the portion of fl + (resp. 
Q h ) consisting of the non-porous rocks. We suppose that 



dn+> £ n <9y e (w) 
a^jn^H = r-» 



We define the fluid part of the domain as 

nj ( W ) = fi+' £ u (r+ u y £ u r-J ( w ) u n-j. 

Let / G C (fi + ), # + G L 2 M 3 ) and #~ G L 2 M 3 ) be functions satisfying 

supp (/) C f2 + and / > in f2 + ; supp (g + ) C f2 + ; supp C f2 A . 
We consider in Q'j (u>) the reaction-diffusion problem with first-order reaction 
-DAu £ + v £ • Vu £ + lZu £ = f in Q £ j (ou) , 

u £ = o on r+ur, 

— - = onS+USfUA.y, 



(5) 



where w e is the concentration of the pollutant, D = D mo i is the molecular diffusion 
coefficient, TZ is a nonnegative reaction coefficient, n is the unit outer normal and v £ is 
the velocity of the fluid, which is the solution of the Darcy-Stokes problems 



/i+(X+)-% M -Vp M = g + infi^, 

/' (K ) ' .<! Vp .,/ 

div (u e>d ) 
-fie 2 Av £yS + Vp e , s 



5 in^>, 
in fi+' e U ft" J, 
on dVL + f ' £ U U r, 



(6) 



in Y £ (u) , 
div (v e>s ) = in Y £ (u) , 
f £jS = on A e (u) , 



with the following interface conditions 
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K,«) 3 = ( v ^) 3 on r J 6 ( w ) u r M ( w ) » 

|*s=0,-h = - Pe,s Oil rj (w) U 1\ £ (w) , 



dx 3 



2 d ( V e,s) T _ sr, + ,-l/2^. ^ j-i-f- ."\ ( 7 ) 



^ 2 ^T^ = "7 (J£T v > e ,) T onlt( W 

^ g = l( K e) (Ve,s) T OnT h£ (uj) 



■,/' 



where: 

• t> £j d and p £j rf are respectively Darcy's velocity and pressure in fit ' e and f2^' 

• ii^ and lf~ are the absolute permeability matrices in f2j' £ and f2^'J respectively, 

• t> £jS and p £jS are respectively the velocity and the pressure of the Stokes flow in the 
fissures, 

• /i + (resp. fx) is the viscosity coefficient in f2^' e (resp. in Oj^'f, 1^ (w)), 

• (^e,s) T = > 0O 2 ) is tne tangential velocity 

We suppose that the 3x3 matrices Kf and K~ are defined through the Z-periodic 
construction: Kf = K + (x/e) and K~ = K~ (x/e), where K + and K~ are bounded, 
symmetric and positive definite, and that /i 1 * 1 and \i are positive constants. In the above 
interface conditions, ©i means the continuity of the mass flux through the interfaces 
Tq £ (cj) and T^ £ (cu), j7J) 2 represents the continuity of the normal stress through the corre- 
sponding interface, and the two last equalities of (J7J) represent the Beavers- Joseph-Saffman 
conditions on the tangential stress, with some nonnegative slippage coefficient 7 (see [3], 
[6] and [33]). 

We first describe the asymptotic behaviour of the solution v £ of ©-© using T- 
convergence methods (see jl] and [8], for the definition and the properties of this varia- 
tional convergence). We prove that the asymptotic velocities v^ d (in fl + ) and v^ d (in fl^) 
and the asymptotic pressures pQ (in f2 + ) and p$ (in f2^) are linked through the Darcy 
laws ( 1251) in Q + and respectively. We also describe the asymptotic behaviour of the 
velocity v £<s (see Corollary [10]). We then describe the asymptotic behaviour of the solution 
u £ of (jSJ) using the energy method. We prove that the flux of pollutant through Tq is 
given through fl39|) ^. while the flux through is given through fl39]U . 

Homogenization theory introduced in the few past decades (see for instance [7] and 
[T2] and the references therein) gives powerful tools for the description of equivalent 
media which are heterogeneous at a microscopic level. The homogenization of transport 
problems of chemical products through porous media has been studied by many authors 
(see for example [11] and [13]). A model of random fissures has already been studied in 
[9] for a problem of radiophysics posed in U (Ujg^Qj (w)) U IR^, where Q £ j (u) is the 
j-th fissure and R^ h = {x G ffi 2 | x 2 < —h}, = {1 G I 2 | 12 > 0}. We here adopt the 
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shape of the fissures used in [9], that we extend to a 3D case, in order to model in a quite 
realistic way the constitution of the soil. 

The paper is organized as follows. In the following section, we introduce and study 
the appropriate local coordinates inside the fissures. In the third section, we study the 
convergence of the velocities. The fourth part is concerned with the asymptotic analysis 
of the contaminant transport problem. We first deal with the case where TZ = 0. We build 
the solutions of local problems in the neighborhood of the fissures in order to pass to the 
limit in the original problem and study their asymptotic properties. The fluxes across Tq 
and T h given in this section are those obtained in (139])^. respectively, with 1Z = 0. In the 
last part of this section, we also consider the case of a dispersive contaminant (TZ > 0) 
with random dispersion in the fissures (see Remark [T71) . We finally give the asymptotic 
behavior of the fully reaction-diffusion problem. We here again introduce the solutions of 
local problems. 



2 Local coordinates in the fissures 

In the fissure Y £ti j (co), G I £ (cj), we define, for a fixed event uj for which the conditions 
([I]) and ([3]) are satisfied, £ x = x± —ie, £ 2 = X2—je, z = —x^, and introduce the curvilinear 
coordinates t G (0, h) and 2/1,2/2 G (—e/2,e/2). Thus doing, the lateral boundary of the 
fissure coincides with the planes 2/1,2/2 = ±e/2. These coordinates are described through 



(£i> £2, ^> 2/1,2/2,*) 

< $2, £ (61, £ 2 , ^,2/1,2/2,0 

$V (61, £2, 2, 2/i> 2/2,0 



= ^-at(e- z) (| + 
= £2 - 4 ( £ ~^) (- + 2/2) 



- a f " (e- e z) (~ - 2/1) = 0, 

- a J ( £ ' 6z ) (2 ~ = °> (8) 
= 0. 



Defining ( 1 = Ci/ £ , C2 = £2/^ an d t = e 9 t, the orthogonality conditions of the 
coordinates given through 

d$i, e d$ 3 , e | d$i, e d$3, e 
and the condition ip £ (0,0,e~ 9 t) = r imply the following Cauchy system 
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d( 2 



which has to be solved. 



-2(1-0) 



a t (i> e ) ~ a i M 

(at (*.))' C '-° r(W 



r 2(l-0) 



- (< We)) 



/ C 2 - s T (VO 



(9) 



r, 



Lemma 1 1. TTie system has a unique solution ip £ (Ci? C25 = r + (Ci, r ) + 
^ e (C 2 , r), u*to ^ (0, r) = V £ 2 (0, r) = 0. 

For every k G N*, every £ 1; C 2 G [ — ^, anc ^ every r£l, one /ias, w/ien e zs c/ose 
to 

V B (Ci,C 2 ,*) = r + O (e 2 ^) , 



<9V e C? <9 2 ?/> £ 



O (e^-Q) 



a = 1,2. 



Proof. 1. Observe that a~[ (ip £ )—al (ip B ) = q {ip £ + otj (w) , w) > Ci > and (a^ (VO) 
C3, thanks to (|IJ). One deduces that, for every I; 6 N* and every Ci, £ 2 G [—k,k], the 
functions 



< 



(Ci,VO 

(C 2 ,V0 



' Ci - 4 WO 



' Cl - «i We) 



! > 



af ty e ) - a; (^J ^ vre " a+ (^) - a" (^) ' 

/ — / / \\' C 2 -<WJ /+ , nV C 2 -Oj wj 



4 (VO - S (V>J ' £ S We) - Oj We) 

are bounded and Lipschitz continuous with respect to (( 1 ,ip e ) and (C^VO respectively. 
The Cauchy problem (J9j) thus has a unique solution ip £ . The condition -?/> e (0, 0, r) = r 
implies V £ (Ci, C 2 , *) = r + ^ (Ci, r) + ^ (C 2 , r), with ^ (0, r) = ^ (0, r) = 0. 
2. Thanks to the preceding point and to the hypothesis (JTj), the quantities 
^ and ^ ^ r are equal to some O (e 2 ^ - ^), a = 1, 2. ■ 

Lemma 2 TTie system (0|) /ias a unique solution (£ x (yi,w 2 , t) , £ 2 (2/1,1/2, , z (2/1, 2/ 2) *)) ■ 
Proof. The Jacobian A of the matrix ^^af^/^-A 3 '^ can be computed as 



A = 1 + 



( ~\i Ci — a i 1 +%/ Ci - ° 



9C 2 



a„- — a« 



aj - 
1 C 2 - a . 



Thanks to Lemma [TJ one has A = 1 + O [e 2 ^ 1 e n, when e is close to 0. One deduces, 
using the implicit function theorem, that the system (JS]) has a unique solution (£i,£ 2 ,z). 



Thanks to the implicit function theorem, one has 



f Ml 

Mi 
Ml 



A v 



2 / + 



(10) 



where 



( a i ) 

(«7) 



;C 2 



4- 



(4) 
(4) 



> Ci - « t 

/ c 2 - a ,- 



One also has 



r % 

dyi 
Ml 

dt 

Let (</ e>c 



2(a+-ar)^ £ 



A 



A 9C 2 

f 1-^.9 \ 

aw 

t 



V 



ac 2 



,9^ 



A 



'a,/3=l,2,3 



o 

dr £ ' v 



dz 
dyi 

dz 
dy 2 

dz 
dt 



A 



a* ) , 



A 
1 



d( 2 
(l 

\ 



(at - a, ) 



(11) 



~<9C 2 



J 



be the metric tensor associated to the local basis defined through 



the vectors fllOp and (jlip . One has the following result. 

Lemma 3 1. The metric tensor {g e>a p) a a =123 satisfies a symmetry property and the 
following behaviour 

9e,n = (ft) 2 (e-°t) + O (e^-V) , ^ £ , 22 = fe) 2 (e^t) + O (s^) , 
9e,i2 = OW-O), ^ >23 = 0(^-*)), 

^,13 = 0(5 2 ^)), ^33 = l + 0(e^), 

where = g (£ _6, t + ctj (u) , co), Vi G Z. 

TTie contravariant components {9sj a p = i2 3 °/ fe,^) a /3=i2 3 satisfy a symmetry 



8 



property and 



9l l 
9? 



O (e 2 ^) 
1 

fe) 2 (e-H) 



— + O (e^-V) , g 



23 

e 



9? 



+ 0(e^~V), g 



33 

e 



(e^-V) , 

1 + O (e 2 ^) . 



Proof. Observing that 



O (e 2 - 3e ) 



det = (g,) 2 fe) 2 (^0 + O (e 2{1 - e) ) 



these formulas are direct consequences of (TTOjl - flTTj) . ■ 

One deduces from the preceding results that the gradient of a function u expressed in 
the local coordinates (j/i,j/2,t) inside the fissure Y e jj (u) is of the form 

1 du 1 du du 



v» = (« + (^))(I^,l 



qj dy 2 ' dtj ' 

for some non-diagonal matrix O (e 2 ^ -9 ^). 

Hence, using the formula Am = "Tpj^a (^\/~\9\9 al3 d pu^j , for a, (3 = 1,2,3, with |g| 
| det (p a/ g)|, one has 



Am 



1 d 2 u 



+ 



1 d 2 u 



+ 



1 d 



(ufdiyi) 2 {q 3 Yd{y 2 ) 2 q^dt 



3 Study of the fluid flow 

3.1 Existence of a weak solution and a priori estimates 

We define the functional space 



V 



veL 2 (Q* f (u) ; M 3 ) | div (v) = in Q.) (u) , v |y eH e H 1 (Y £ (u) ; R 3 ) 



v = on A £ (u) , w • n = on <9^' £ U U T 

V £ is a complete space when endowed with the norm 



II 2 n 



+ \\v\ 



f 



(w);K 9 ) 



1/2 



9 



Multiplying fl6]) by $ e V £ , using Green's formula and the conditions ([7]), we obtain 
the following variational formulation 

H+ (K+)' 1 v £>d ■ <f>dx + fx' / (K-y 1 v e4 ■ $dx + fie 2 Vv £ , s -V$dx 

Jttp £ ^ n h,'f JY e (u) 

+7 / (^ e + )" 1/2 («.,.) T • ($) T dx* + 7 / (^ £ -)~ 1/2 (««,.) T • (*) T dx' (12) 



g + ■ &dx + g ■ <f>dx 



where x' = (xi,x 2 )- Using standard arguments, we immediately deduce from this vari- 
ational formulation that the system flSJ)-([7|) has a unique weak solution (v £ ,p £ ) G V £ x 
L 2 (tt £ (w)) /R. 

Lemma 4 There exists a non-random constant C independent of e such that 

2 9 / I 2 / ■ 1 2 

|t> £ | <ix < £7e / |Vf e | dx, / |t> £ | dx < C, 

>Y e (w) Jy £ (lu) Jn £ f (uj) 

\v £ \ 2 dx'+[ \v £ \ 2 dx' < Ce 2 , e 2 [ \Vv £ \ 2 dx < C. 

■o, e ^r- £ yy e ( u ) 

Proof. Define the normalized fissure 

7 , n _ / (21,22,^3) e R 2 I or (-e - "^) < 21 < a+ (-e~ e x 3 ) , 1 , , 

Z " U ' ' l »- (-<r fl * 3 ) < z 2 < a+ (-rt 3 ) , x 3 G (-/», 0) r ' 1 ' ° ' 



a 



For every f 6 C 1 (-^ £ ,y (w)) such that j/i = on the lateral boundary of Z £)i j (a;), one 



has inside Z £>i j (u) 



^ (s ' z2 ' xs))2= {lM dzi ) - q ^ £ " x3) I a ~ (H) dzi - 

Thanks to the hypothesis (CQ), one has the following Poincare estimate 



ip dz\dz 2 <C |V^| dz\dz 2 . 

Defining z\ = (x% — si) fe and z 2 = (x 2 — ej) fe, one gets the first estimate. 

We replace $ by v £ in the variational formulation (|12p . Since the matrices Kf and 
K~ are bounded, symmetric and positive definite, we deduce, using the first part of this 
Lemma, that 

/ \v e4 \ 2 dx + I \v E ,d\ 2 dx + [is 2 I \Vv £ , s \ 2 dx + 7 / \(v £tS ) r \ 2 dx' 
+7/ \{v £ ,s) T \ 2 dx' < C \ / g + -v £4 dx+ / g~ ■ v £>d dx 



07 
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hence, using Cauchy-Schwarz inequality 

1/2 



\v e \ dx < C 



and 



\v £ \ dx I + C 



dx Y iC {L, 



1/2 



\v e \ dx 



(w) 



/ |Vf £ | dx < C I / |f e | z (ix 



1/2 



Using a trace theorem in the normalized fissure Z £<i j (u), we have 



ip 2 dz 1 dz 2 < CI | Vip\ 2 dz\dz2i 

which ends the proof. ■ 

We now deal with the pressure. We define the extension p £ ^ of the pressure p £ ^ in the 
solid parts of the porous media 

{Pe,d in fi+' £ U tt~ ,e f , 

t4iT / Pe,d(a £ z + a £ l)dz Vz E eZ 2 + le C tt+ U tt h , I e Z 3 . 

We define the zero mean- value pressures 

1 



Pe,d 



PI 
Pi 

We = Pe,s 



Pe,ddx 



Pe,d - 77—7 / p e ,dda; 



in VL+ 
in f2 



/l ! 



Lemma 5 There exists a non-random constant C independent of e such that 

[ \Vpt\ 2 dx < C, [ \Vp £ \ 2 dx < C, [ {Wfdx < C, 
Jn+ Jn~ Jy £ (lu) 

[ {ptfdx < C, j {p £ fdx < C. 

h 

Proof. We multiply (JE])i,2 by Vp e and obtain 



£ ,ddx, 



/ \Vp £ , d \ 2 dx = - g + ■ Vp £ , d dx + / (Kf) 1 v £ , d -Vp 
Jnj' e Jn+' E Jnp £ 

I \Vp e4 \ 2 dx = - g~ ■ Vp e4 dx + n~ / {K-y 1 v e4 ■ Vp e4 dx. 

Jnr'j Jo.7' e , JciZ'i 
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The boundary conditions (JEJU, the smoothness of and Lemma |4] imply 

\Vpt\ 2 dx<C; [ \Vp~\ 2 dx<C. 



Using Poincare-Wirtinger' inequality, we have 



{ptfdx < C(fl+) / \Vp £ \ 2 dx < C, 

Jn+ 

(p-fdx < C(n h ) / \Wp £ \ 2 dx < C. 



m- Jn h 



In order to get estimates on p £ , we consider the problem 

div 2 ($°) = %(ez U £z 2 ,x 3 ) in Z £i!J - (w) , 
$ £ ° = on dZ E ^ (u) , 

where the fissure Z £ ^j (u) is defined in (Fl3l) . This problem has a unique solution $° in 
the space {(— A) -1 Vw \ w G L 2 (Z £ ^j (w))} (see [17] for example), such that 

f |V$°| 2 dzdx 3 < C j (p- £ ) 2 dzdx 3 , (14) 

where, due to the hypothesis ((!]), the constant C is non-random and independent of e. 
Let us define 

/ s ($°) 1 ((xi - is) /e, {x 2 - je) /e, x 3 ) 
$ £ (x) = £ ($°) 2 ((xi - ie) /e, (x 2 - je) /e, x 3 ) 
\ ($°)s(( x i -ie)/e,{x 2 - je)/e,x 3 ) 

Then 

div ($ e ) (x) = div 2 ($°) (x) = p7 (x) in y £iij (w) , 
$ £ = on 9r e , y (w) 



and, thanks to (EE 



/ |V$ £ | 2 cfe< / |V$°| 2 ^x 3 < C / (p- £ ) 2 dzdx 3 = ^ [ (p- £ (x)) 2 dx. 

(15) 

Multiplying ©5 by $ £ in F e (w), we get 

2 / V7„. V7/F> J™ / /'^-/'™\\ 2 



/ie 2 / Wv £ ■ V<S> £ dx - / (p7(z)) cb = 0. 



Using the inequality flTojl . we have 

(Pe ( a; )) 2 ^ < C ( / (pr(^)) 2 ^) ( £ 2 / \Wv £ \ 2 dx 
Thanks to Lemma HJ we derive the last estimate in Y £ (u). ■ 
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3.2 Convergence 

Observe the following result. 

Lemma 6 1. For every ip e Cq (Y^), we have 

lim / ip (xi,X2,Xz) dx = h (q 2 ^ / ip (x', 0) dx' , almost surely. 

2. Let {w £ ) £ be a sequence such that sup e j Y , , (u> e ) 2 <ix < +oo. T/iere exists a subse- 
quence, still denoted in the same way, such that 

lim / w £ (pdx = h (q 2 (0)) w (x', 0) (p (x', 0) dx', VV> € Co (M 3 ) , almost surely. 

Proof. 1. Define £ ie (y) = ei — eri (e~~ e t) — yiq-i (£~ e t) and (y) = ej — £Tj (s~ e t) — 
yiqj {s~ e t) where r { (s^ e t) = r [e~ e t + /3 { (u) , u) and q^ {s~ e t) = q (e~ e t + a { (u) , u). 
According to the properties of the above-defined curvilinear coordinates, we have 

Km V f 

(ij)6/e(w) e '* J ( ' 

fe/2 re/2 rh 

= liH1 n 22 / / / ^(^(f).^^).-*)^^)*^)^ 



if (ei, ej, -t) q 2 (e^ 9 t) dt 
lim _ / I / <p (xi,X2, — £ 9 s) dx\dx2 I g 2 (s) ds 



/i (g 2 ) / p (x', 0) da;', 



s 

where we have used the ergodicity result (j2J). 

2. The sequence of measures (V e ) e , with z/ e = lw^dx, 1^ being the characteristic 
function of the set A, thus converges in the weak sense of measures, when e goes to 0, 
to the measure v = h (q 2 ) l s (x') dx'. Using the hypothesis on (w £ ) £ , we deduce that 
the sequence of measures {y £ w £ ) £ has bounded variation. Up to some subsequence, the 
sequence {v £ w £ ) £ thus converges to some Xo i n the weak sense of measures. For every 
if G Co (K 3 ), one has, thanks to Fenchel's inequality 



/ {w £ f du £ >2 w £ ipdu £ - / ip 2 dv £ . 

JR3 Jm.3 J R 3 

Passing to the limit, we get 

+oo > liminf / (w £ ) 2 dv £ > 2 (x , <p) — I f 2 (x', 0) dv. 
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Thus 



sup I (xo, <p)\<peC (R 3 ) , J^y (x', 0) dv < 1 1 < +00. 
Using Riesz' representation theorem, we can identify x with wu, for some w £ L 2 (R 2 ) 



Remark 7 i. From Lemma 0, nsmg £/ie ak>i>e result, we deduce that, up to some 
subsequence 

lim / v £ -<$>dx = h (q 2 (0)) / v f ■ ®dx' , V$ G C (E; R 3 ) , 

almost surely, and from Lemma\f^ we deduce the existence ofiiQ G L 2 (E) snc/i i/iai, 
to some subsequence 

lim / ~p~ £ tpdx = / -Koipdu = h (g 2 ) / 7r (x') y2 (x') dx' , Vy? G C (E) , 
e ^°iy £ (aj) Js is 

almost surely, where x' = (xi,^)- 

In i/ie resi 0/ the paper, we will no more indicate this almost surely convergence 
where there is no doubt. 

2. From Lemma\J^and the above computations, we deduce the existence ofvo G L 2 (Q; R 3 ) 
such that, up to some subsequence 

v £ \ n +,e ->• v | n +=: v+ d w-L 2 (fi+;R 3 ), 



v * in 



->• v \ n -=:v od w-L 2 (£l h ; R 3 ) 

lim / v £ ■ §dx = h (q 2 ) / v ■ §dx' V$ G C (S 

£ ^ JyJlo) Jt, 



5. VFe set t>o |e=: t>o,/- For ever?/ <p E Cq (£), one /ias 

= lim > / div (n e ) o?<ir = —lim > / t> £ • Vwdx 

(i,i)6/ e (u)' / y <W^ (j,j)G/ e (aj)^ y <W W 

= — /i (g 2 ) v Q f Wipdx' 

r ' 

= h(q 2 ) / div (vqj) pdx', 



thanks to the estimates of Lemma^ Thus div (i>o,/) = in E. 

^. It zs easzfo/ seen that div (n^) = in fi + and div (n^ d ) — in f^~. For ever?/ 
G C 1 one /ias 

= lim / div (v £ ) tpdx = —lim / v e ■ V^pdx + lim / (v £ )n ipdx'. 
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Observe that 



lim 

£->0 



{v £ ) 3 <pdx' = (q 2 (0)) / {v 0J ) 3 ipdx' 



whence {vq^) |sx{o}= (q 2 (0)) (^o,/) 3 - I n a similar way, but working in Q h j, instead 
of tip 6 , we have (v^ d ) 3 \zx{-h}= (<? 2 (°)> ( v o,f) 3 - 

5. From Lemmas{5\ and® we get, up to some subsequence 

^+ - «+ s-L 2 (fi+), 

%tpdx = h (q 2 (0)) / TT (pdx' \/(p £ C (E) . 



PI -> 
Pi ^ n Po 



lim 

e-i-0 



(16) 



We now set 



Vr 



v £ L 2 (n+ U R 3 ) | div (v) = in ft+ U 

« • n = o on r + u r~, v 3 |ex{o}= (q 2 (o)) (v,/) 3 = ^ 3 |s X {-m 



Every function t> £ Vo can be extended in a function of L 2 (y e (u) ; IR 3 ) independent 
of x 3 in F e (tu). 

We define the appropriate notion of convergence for the problem (JH]). 

Definition 8 A sequence {V e ) £ , with V e £ V e for every e, r -converges to some V £ V 
if 



f e->0 



e in 



V 



w-L 2 (fi + ;M 3 ) : 
w-L 2 (n h ;M. 3 ) 



lim / V £ ■ $dx = h (q 2 (0)) / V f ■ ®dx' V$ £ C (E; M 3 ) 
Vf := V |s. 

We define the functional F £ on L 2 (tu) ; IR 3 ) through 

fi + {Kfy 1 v ■ vdx + fi (K^y 1 v ■ vdx + fie 2 |V 

Jn+' e J Q hj Jy e (uj) 

+7 / (K+)' l/2 v T ■ v T dx' + 7 I (K-y 1/2 v T ■ v T dx' 



FAv) 



v\ 2 dx 



-<x> 



if v £ V £ , 
otherwise. 



In order to describe the asymptotic behaviour of this functional, we consider the Z- 
periodic solution $^ of the local Darcy systems 

(i^r 1 ^ - Vtt± = e k inZ\ k = 1,2,3, 
div(<l>±) =0 inZ 1 , 
• n = on S 1 , 
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where (e±, e^) is the canonical basis of M 3 . We then consider the Stokes system 



-A Vk + V^ k = e k mZ' = (-l/2,l/2)\k = l,2, 
div (r] k ) = in Z', 

Vk 



(17) 



on dZ' 



where (ei, e-i) is the canonical basis of IR 2 , and the local scalar problem 

\2 



1 in Z' = (-1/2, 1/2)' 
on dZ'. 



We define the 3x3 matrices K + , K , the 2x2 matrix Kt and the constant ko through 



K 



ml 

(K f ) ml 
k 



i dz 

i ( $ m)/ dZ 
/ (Vm)ldz 

| Vrj \ 2 dz = t] Q dz. 

1 JZ' 



m, I = 


1,2,3, 


m, I = 


1,2,3, 


m, I = 


1,2, 



(19) 



One can prove that the matrices K and Kf are symmetric and positive definite (see 
16]). Our main result of this part reads as follows. 



Theorem 9 Suppose that r is also a stationary random process. The sequence (F £ ) 
Y -converges in the topology tq to the junctional Fq defined through 



Fa(v) 



K+ ) v d ■ v d dx + V 
h (g 2 > 



K~ 



v d ■ v d dx 



(?) 



+ <9 2 > 7 / ((A- + r 1/2 + (K-y in ) (v,) T ■ (v,) T dx' 



+ 00 



where K*± = ffc-ff /A- f S K± (*i> z ^ °) dz x dz 2 . 



((v f ) 3 ) 2 dx' 



ifve V , 
otherwise, 



/<o-(0)> J(a-(0)) 

Proof. For the definition and the properties of the T-convergence, we refer to [I] and [8]. 



Upper r-limit. Choose any smooth v G C 1 ^f2 + U Q 
■ Vj, v |e=: Vf and build in Y Eti j (to 

((K f )- l V f ) k (lE,]E^) f° 



n V . We define v 



«/)s 



h(q) 

(vf) 3 (ie,je,0) 

hkn 



Ve,k,ij ( X ) dx ^ 



(20) 



VefiAj ( X ) dx 3 
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where 



with 



V e ,o,ij( x ) = r ]o (z(x 1 ,x 2 ,x 3 )) , 



z(x 1 ,x 2 ,x 3 ) 



( xi-ie-e(a i (-£ e x 3 ) + a+ (-£ e x 3 )) /2 \ 

/ 



eqi (-e- e x 3 ) 
x 2 -je-e (aj (-e" e x 3 ) + at (-e^x 3 )) /2 



?i ( s ) = 9 ( s + a i ( w ) , w ) > 



eg,- (-£" e £ 3 ) 



(21) 



r] k being the solution of (ITTj) . ??o the solution of ( TT8|) and Ify and fco being defined in (fl9l) . 
We then define the test-function w° through 



-1 



-1 



7 ,o _ 



in fi+' e , 



in n h j, 

in K M • 



(22) 



One deduces from this construction that v° e \ n +e L 2 (f^' e ; M 3 ), v° | n -G L 2 R 3 ), 

u° |y e ( w )G H 1 (y e (w) ;M 3 ), f ° = on A e (w), div (v°) = in fi+' £ U fi^ e , w £ °-n = 0on 
<9f2j (u) and 



div (t; £ ) 



h(q) 



f° 

/ div 2 (^.Mi) ( x ) dxs = °> in Y *M 
J-h 



Therefore f° G V e . Moreover w° is independent of x 3 in each fissure Y E} ij (u>). Using 
the ergodic result ([2} and making some computations, one easily proves that (i>°) E To- 
converges to t>. 

We compute the limit \im E ^, F e (i>°). We have 



lim/xe 2 / |Vt>°| 2 Gfe = lim/ie 2 / |V(u°) | 2 cte 

+limu£ 2 £ / |V(t£) s | 2 dx. 



Using the expression (1201) 1 . we have 

E / y |V(^)J 2 rfx = lim / . 2 £ * 

/"° 

((^/) 1 «/) fc (ie, je, 0) / qi(-e~ 6 x 3 )V z r) k {z{xi,X2,x 3 ))dx 3 



( 



Ye,ij{u) 



\ 



x ((K f ) 1 v f ) l (ie,je,0) / ft (-£ e x 3 ) (2 (a?i, x 2 , x 3 )) cb 3 



/ 
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We introduce the change of variables (zi,z 2 ) = z (x%, x 2 , xs), where z (xi, x 2 , £3) has 
been defined in (l2T]) i . and get 



lim/ie 2 / |V (u^^da; 

= (ilim v 4/V f (TO- 1 »,) ? (feo>.o)((^)-'» / |,(fe, J >,o) N 

e_>0 (iti fe (w) (Q) J-hJz' \ x ^zVk (z) ■ V z m (z) q { (s °x 3 ) qj (-£ y x 3 ) d^dar 3 / 



Using the ergodicity result (jSJ) and the definition (JT9J) 1 of if/, the above limit is equal 

to 

(?) 

because Vf is independent of x 3 in Y h . Using a similar argument, we have 

(? 2 > (I/? 2 ) 



On the other hand, observe that, for every ■?/> 6 (fi; M?) we have 
lim / v^-ipdx = / ((-^ +N ) uj| / ®f (z) dz ■ i/;dx 



vt ■ ipdx, 



d 

4(0) /-at(O) 



lim/ AT+ (a) • t/j (x) dx = lim > e 2 / / if + (?/', 0) dy' ■ i[> (ie, je, 0) 

(a+(0)) ^(o+CO)) \ 

fT + (21, 22 j 0) dz\dz2 ] ■ V* 0) da/. 



S \J(a-(0)) ^(a-(0)> 



We thus obtain 



lim 



(p+f (Kty 1 v° £ ■ v° £ dx + f,- [ (K-y 1 v° e ■ v° £ dx \ 
Jn+ Jn~ 



+7/ (K+r 1/2 (v° E ) T -(v £ ) T dx' + 1 (K £ r 1/2 (v £ ) T .(v £ ) T dx' 



= fi+ / ttf ■ u+dr + fi f(R ^ v d ■ v d dx 

+ (? 2 > 7 / ((^* + )- 1/2 + (K*-y 1/2 ) (v f ) T ■ (v f ) T dx', 
whence lim^o F £ (v £ ) = F (v). 
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For every v G V , there exists a sequence {v m ) m C C 1 U i\ ; IR 3 j fl V such that 

v m -» u, s-L 2 (n + U 07"; M 3 ) . (23) 

m— >+oo 

Building the sequence ((v m )° £ ) £ associated to v m through fl22l . the sequence ((f m )°) 
ro-converges to v m , and using the above computations for smooth functions, we have: 
lim e _» -F e = F Q (v m ). Hence lim™-** lim^o F e ((f m )°) = -F (u). Using the diag- 

onalization argument of jU Corollary 1.18], there exists a sequence (i>°) e , f° = (t> m (e)) £ 
(m (e) — 7- £ _+o +°°)) sucn that (f £ ) e r -converges to v, and lim sup^o F £ (t>°) < F (v). 

Lower T-limit. Let (i> £ ) £ be a sequence such that t> £ G V e for every e, and (i> £ ) £ 
ro-converges to v. We write the subdifferential inequality 

/i£ 2 / \Vv l £ \ 2 dx>^6 2 [ \V{v m f £ \ 2 dx + 2fie 2 [ V (v m )° • ( V^ 1 -V(v m )°)dx, 

(24) 

where ((f m ) £ ) £ is the sequence associated to v m through (1221) and where the sequence 
(v m ) m satisfies the conditions ff23l . Observe that 

/ V (v m )° e ■ (Vvl - V (v m )° e ) dx = - A (v m )° £ ■ (yl - (v m )° £ ) dx, 

because v\ - (v m )° E = on dY £ (w) \ (r+ U T h £ ) ( w ), and 2{g£ U =0 = | X3= _ fc = 0. 

Then, using the ergodicity result (j2J), we have 



lim2/^ 2 / V (v m )° • (V^ 1 - V (v m )° £ ) dx 



whence 



_ h(q)(y/q) IJ^A zVk (z) {(Kj)- 1 (v m ) T ) k ■ (v - v m ) T (x') dzdx' 
-fi h ^ q ^yi q - I I A z i] (z) (v m ) 3 (v - v m ) 3 (x') dzdx', 



Y.J Z' 



lim lim2/i£ 2 / V (v m )° £ ■ Nvl - V (v m )°) dx = 0. 

Recalling the inequality (1241 and the computations built in the above case of smooth 
functions, we have 



liminfue 2 / \\7vlf dx 

^° J Y e{u) 



{q) Jy, fc o Jt, 



Thus, computing in an easy way the liminf of the other terms in F £ (u £ ), we obtain: 
liminf.F,^ 1 ) >F (v). m 

Let us write the problem associated to the limit functional F Q . 
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Corollary 10 The solution (v £ ,p £ ) of the problem verifies the following properties: 

• (v e ) e r -converges to v G V , and set v \ a += v+ d , v \ n -= v~ d , v | E = v 0J . 

• (pf) £ converges to Pq , s-L 2 (Q + ), (p £ ) £ converges to p$ , s-L 2 (f^) and [Tfyh holds 
true. 

• v od> v od> v o,f; Po ; Po anc ^ are solutions of the following problems: 
i) in the regions Q + and Q^, one has the Darcy laws 

fi + (K + y 1 v+ d -Vp+ = g+ mn+, 
div (vQ d ) =0 in Q + , 

^ {k-^i v- d -Vpo = g- inn-, (25) 

div (v) = mn h> 

i v o,d) 3 lsx{o} = (q 2 ) K/) 3 on E x {0} , 
fc) 3 lsx{-M = (g 2 )K/) 3 Ex{-/i}, 

m ) on S, i/ie velocity (fo,/) 3 given through 

(vo,f) 3 &) = {pt W, 0) - Po (x>, -h)) ^ (g2(0) ^° (1/g2(0)) (26) 



and the tangential velocity (voj) satisfies the modified Darcy law 
(K f rHvoj) T + l((K* + r 1/2 + (Kn- 1/2 )(voA + Vn 



in E, 



div (v ,f) T = in E, 
( v o,/) r " n — on <9E. 



Proof. Thanks to the properties of the r-convergence, the sequence (v e ) £ ro-converges 
to Vq G V and lim^o F £ (v e ) = F (vq), where v is the minimizer of the problem 

inf I Fq (v) — 2 / g + ■ vdx — 2 / g~ ■ vdx 1 . 



For every V G Vo, we have the following identity 

H+ I "* v 04 ■ Vdx + fi~ I (R-y 1 v 04 ■ Vdx 

+p;h^[ (Kfr'M^Vrdx' + fih^p-il/q 2 ) [ (v 0J ) 3 V 3 dx' 
{q} Jt, fc o Jy, 

+ {q 2 ) 1 j^ ((^+)" 1/2 + (K*-y 1/2 ) (v 0J ) r ■ V T dx' 

g + ■ Vdx + I g~ ■ Vdx. 
<>| •/<>, 
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We infer the existence of a pressure p£ (resp. p Q , 7r ) in Q + (resp. Q h , S) such that 

/ Vp+ ■ + / Vpo - Vdx + h [ Vtt ■ K-dx' + \ih ^ ^ ^ Q ^ [ (v 0J ) 3 V 3 dx' = 0, 
Jn+ Jn~ Jt, Jt, 

which implies, because V (resp. V T ) is divergence-free in Q + U (resp. X) 

<q 2 }(l/q 2 } \ , f 

(po 0) - p (^'j _ ^)) + M ; ( v oj) 3 V 3 dx' + h 7r n • V T da. 

's V fc o / Jas 

This gives the result. ■ 

4 Study of the transport problem 

Let us now consider the transport problem (jSJ). In this section, we will describe the 
asymptotic behaviour of the solution u £ of (jSJ), when e goes to 0, distinguishing between 
the cases 1Z = and 7£ 7^ 0. 

4.1 Existence of a weak solution and a priori estimates 

We define the space 

#r +ur - (0/ (<*>)) = {HGiJ 1 (fij (w)) I w = on r + U r~} . 

Lemma 11 i. T/ie problem (TJP /ias a unique weak solution u e G H^ +ur _ (fE (a;)) 
which is nonnegative in Q e j- (00). 

2. There exists a non-random constant C which is independent of e such that 

{u £ fdx <C ; I \Vu e \ 2 dx < C. 

3. There exists a linear and bounded extension operator P £ : H 1 (fH (a;)) — >■ H 1 (Q) 
and two non-random positive constants C\ and C2 such that 

P £ u £ = u e in Q £ j- (co) , 

\u £ \ 2 dx, 

I Vu £ \ 2 dx. 









< 


Cl 


1 








< 




1 









21 



Proof. 1. Using the standard variational methods, one proves that the problem ([5]) has a 
unique solution u £ G H^+ ur - (Oj (w)). Multiplying ([5]) par (u £ )~ = min (0,w e ) and using 



Green's formula, one has 



D |V(u £ ) | 2 rfx+ / (f £ -V(w £ ) ) (u £ ) dx + Tl ((u E ) ) 2 dx 

= [ f(u £ ydx<o, 



because / is nonnegative. Because div (v £ ) = in Q e j (u) and v £ ■ n = on <9£7^ (co), one 
has 

/ (f e • Vw e ) M e C?X = — I (v £ ■ Vw e ) M e C?X = 0. 

Jn e f (u) Jn £ f (uj) 

We deduce that J^e/^ | V (w e ) _ | 2 dx + 7?. ((w e )~) 2 <ix < 0, thus (u £ )~ = 0, in 

f2j (u;), hence u £ is nonnegative in [yS). 
2. As already observed, we have 

D l | Vu £ \ 2 dx + 1Z I (u £ ) 2 dx = / fu £ dx. 

For 1Z 7^ 0, one deduces from this equality, using Cauchy-Schwarz' inequality, that 
/nj(«) ( M -) 2 ^ < C and |Vu £ | 2 da; < C. 

In the case 1Z = 0, one can prove, using [21 Lemma 3.4], that there exists a non-random 
constant C independent of e such that 

u F ) 2 dx < C / |VuJ 2 dx. 



e 

Thus, using the above equality, we get the desired estimates. 
3. This is a particular case of the result given in [TJ. ■ 

We will still denote by u £ its extension P e u £ to the whole Q. From Lemmas dU and 
IHl we deduce the existence of Uq G and Uq G (^), such that, up to some 

subsequence 

u £ | n+ -»> m+ w-#* + , 



lim / u e (fdx = h (q 2 (0)) / u+ (V, 0) 9? (>') dx' Vy> G C (£) 



'y,(w) 

We intend to describe the problems satisfied by Uq and Mg , in their respective domains. 
We now define the notion of convergence associated to sequences satisfying the above 
convergences. 

Definition 12 A sequence {U £ ) £ , with U £ G ifp +ur _ (f^ (a;)) for every e, Ti-converges to 
U, with U \n+=: U + G {Vt + ) and U \ n -=: U~ G H^_ if the convergences (Ffy 

are satisfied, replacing u £ by U £ . 
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4.2 The asymptotic behaviour in the case 1Z = 

In this subsection, we deal with the case 7Z = 0. Using the boundary conditions ©2,3, 
we consider the variational formulation of the problem (J5]) 

Vu G ifp+ ur _ (cu)) : D / Vu £ ■ Vudx + (v £ ■ Vu e ) udx = / fudx, 



where t> £ is the velocity of the fluid flow, that is the solution of (E])-([7j). We consider the 
Z-periodic solution bj of the cell problem 



Abj = inZ 1 j = 1,2,3, 
(V&j + ej) ■ n = on 5 



(2f 



Ac m = in Z' m = 1, 2, 



and the problem 

{ (Vc m + e m ) • n = on 9.Z 7 . 

We define the tensors D and £)* through 

- / r db- 

Dij = D\\Z l \5 i:j + Ig^dz 

Let x+ (resp. xi) be the characteristic function of f^' 6 (u>) (resp. fi^'J (a;)). We have 
the following result. 

Lemma 13 One has, up to some subsequence: 

1. x £ +DWu £ ^ £ ^o DS7u^, in L 2 R 3 )-weak, 

2. x £ -DVu £ ^ £ ^o DVu , in L 2 (Q^;R 3 )-weak, 

3. lim £ ^o j Ye{u) DV T u £ ■ <pdx = h(q 2 (0)) / E D*V T < (a/, 0) • <^r' ; G C~(S;M 2 ), 

Proof. 1. Let tp G C c °° 6^ = ebj (x/e). Multiplying © (for ft = 0) by X+&j¥>, we 

get 

/ x+DVu £ ■ (<pVbj + bjV(p) dx+ x e +v £ ■ Vujf^dx = / x e +fbfodx, 
Jn+ Jn+ Jn+ 

from which we deduce that 

limy x £ +DVu e cp ■ Vb)dx = 0. (29) 
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Observe now that, through (1251) 

/ x+DV (vu £ ) ■ (Vfr; + e 3 ) dx = D x% (V6^ + e 5 ) <p ■ Vu £ dx 
Jq+ ' Jn+ 

+D / x £ + (V^ + ej) u e ■ V^dx 

= 0. 

Thus, taking into account (129|) . one has, up to some subsequence 
lim / x+D(p [e 3 - ■ Vw £ ) dx = — MmD / x+ (ytf + e j) u e ■ Vyjcix 

- -»£(i*iW,£*)/.&** 



2. In a similar way than above, we get, for every <p G (Q h ) 

3 



lim / X-Dif (e,- ■ Vu £ ) dx = Dj^ (\ z1 \ &ij 
3. We define the quantity cL, m — 1, 2, through 



9m 



<9x 



—ipdx. 



c e m {x) =e(q (0))c ri 



( x 1 -ie-e(a i (-£-°x 3 ) + a+ (-£- x 3 )) /2 \ 

/ 



x 2 -je-e (aj (-e~ e x 3 ) + aj (-e~ e x 3 )) /2 



Then, using curvilinear coordinates, the ergodic result ([2]) and making some compu- 
tations as before, we get the result. ■ 

Our main result in this subsection reads as follows. 

Theorem 14 The sequence (u e ) , where u £ is the solution of fifty, converges in the topology 
T\ to the solution (u~q,Uq') of the variational formulation 



V (it+ vr) E Hl + (n+) x (n h ) : / DVu+ ■ Vu+dx + / DVu^ ■ Vurdx 

Jn+ Jq- 

+ / i v o,d ■ V«o ) u+dx + / ( v o, d ■ Vu o ) u ~ dx 
Jn+ Jq- 

+h (q 2 (0)> / D*V T u+ ■ W T u + dx' - h (q 2 (0)> / ((v ) T • V T w+) u+dx' (30) 
fu + dx, 



D 



h(l/q 2 )h 
= \Z l 



Po ~ Po k o 
D(q 2 ) (1/q 2 ) n 



dx' 
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where v^ d and v d are the limit velocities appearing in Remark [?| and p$ and p are the 
pressures appearing in CorollaryUU 

Before starting the proof of Theorem [T3J let us introduce the constant "vertical" 
velocity (i> e ,y) 3 in the fissure Y £yi j (a>), G I £ (u), defined as 

( V £,ij)3 = {Pa,e,ij ~Po,e,ij) ^ / 2 \ l\ I q^V ^ 

where Po e>i j = Pq (i£,je, 0), Pa~ ei j = Po (i e i3 e i Po and p^ being the pressures defined 
in Corollary [TOl (compare to ff26l) ). Inside the fissure Y £ ^ (co), for every G I £ (cj), we 
define, for every u G C 2 (Q) satisfying u = on T = T + U T~ 



u £)ij (xi,x 2 , x 3 ) = u + (xi, x 2 , 0) 

, , :i ^(---"n^(--- fl , ' M ' V £> 



o ^ f t{v £i ) 
(u~ (xi,x 2 , -h) - u + (xi,x 2 ,0)) / — ? — ^— — ? — ^— -expl 6 -^-^- ) dt 



+ 







exp — - dt 



- h *(-e-H) qj (-e-H) * \ D 
where g« (s) is defined in ( |2Ti) 9. We finally define the test-function m , £ through 



(32) 



u + (x) in 

Mo, £ (^) = { u £ ,ij (x) in K £)i j (w) , V (i, j) G J e (w) , (33) 
m _ (x) in f2^. 

The properties of this test-function are gathered in the following result. 
Lemma 15 1. One has: 



(-e -0 ^) (« 6 ,v) 3 -7^- ( s ) = in F Mj (id) , 

w £)ii (xi,x 2 ,0) = u + (x 1 ,x 2 ,0), onT^ eij (u) 
u £ ,ij (xi,x 2 , -h) = u~ (xi, x 2 , -h) on T~^ Aj (u) 



(34) 



2. For every e > 0, u Q>£ G H^ +ur _ (fl £ (a;)). 

3. The sequence (wo,e) e T\-converges to u. 

Proof. 1. This is an immediate consequence of the definition ( 1321) of u £ ,ij- 

2. This is an immediate consequence of the construction ( 1331) of wo, £ , m ^ + and in 
and through the "boundary conditions" (1341) 9 3 satisfied by u 0t£ at the ends of the fissure 



25 



3. From this construction, we deduce that (uo )E |n+) e (resp. (uq >£ | n -J ) converges to 



u (resp. u ) in H 1 (fi + )-strong (resp. H 1 (Q h ) -strong). 
Moreover, for every f e C (M 3 ), we define 

cO 

exp 



One has 
limA F 



e^O 



and 



limAg 

£^0 



I ^ 



D 



-hQi(-e- e t)q j (-e-°t) 



exp 



lim > 

£->0 



(i,j)ei e {u>) 



h qi(-e~H) qj (-e-H) 
( f° 1 



exp 



x I if (ei,ej,x 3 ) 

-h 



D 
\ 



X3 qi (-e-H) qj (s-H) 
x exp 



*K«) 3 i (/f 



g 2 (— £ e x 3 ) <ix 3 



lim > 

e->-0 ^— ' 



(ij')e/ e (w) 







ft (~e~**) 



exp 



\ 



9fgi(-e- fl *)^( T e- 9 t) 
x exp 



<? 2 (0 ^ 



Thus 



£^0 



X 



lim / </? 



u (xi,x 2 , —h) -u + (xi,x 2 ,0) 



o 



- h qi(-e-H) qj {-e-H) 



exp 



' X3 qi(-£- e t)qj (-e-H) 

from which we deduce that 



exp 



^ 3 i dtdx = 0, 



£->0 



lim/ ifUo^dx = h (g 2 (0)) u + (x , 0) f (x , 0) dx' 



Y e ( U ) 



Thus the sequence («o,e) s Tx-converges to w. ■ 
Proof of Theorem 1141 Thanks to the boundary conditions ©4 7, we have 



/ (v e ■ Wu £ ) u 0:£ dx = - / (v e • Vw ,£) 
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Using Lemma [T5l and the "compensated compactness" result (see [H]), we immediately 
deduce the following limits 



ton/ (v t .Vu M = / w ,v^v*. 

£ ^°in+ Jn+ 

lim / (t> £ • Vw e ) u 0j£ dx = / (t^ d • Vmq ) u~dx, 



Uo e being independent of e in Q + U fL . We then write 



D / Vu £ ■ Vu , £ dx - I (v £ - Vw , £ ) u £ dx 

= D V T u £ ■ V T u 0)£ dx - / ((v £ ) T ■ V T M , e ) u £ dx 

JYJui) JYJuj) 



+D 



du £ du 0:( 
y s (uj)9x3 dx 3 



■dx 



Y e (w) 



( r ... ) 3 ^^u £ dx. 
1 ox 3 



A direct computation gives 
du £ Ouq 



lim ( D 



Y E {u) ox 3 9x 3 



■dx 



du . 



lim > 



du F du 



v £<s ) 3 -g— ] u £ dx 



du FAA 



c \77 1 \ J Yeij(u) \ 9x 3 9x 3 
( ^ f ( n 5 



D^—^-qiqj - (f £)S ) 3 -^r~u £ qiqj ) dx 



9x s 



lim 

e-s>0 



E 



9x3 
9m 



9m 



du £tij 
dx* 



u F dx 



, ..(a;) ^3 

9m 



^J^- |x3=o <Mj (0)u £ dx' 



E 



(i,i)e/ E (w)" r fc. E .ii (aj) 



9x a 



\X3=— h 



QiQj (—he e ) u £ dx' 



Hence 



lim ( D 



du £ du 0:l 

Y e (fj) 9x 3 9x 3 



dx 



lim 

e-H) 



E , 

(tj)6/«(w) ,/ye >«^ 



9m 



(( u £,*i) 3 - 00 3) -Tr^qiQjUedx 
0X3 

■ D E / r+ M (*»(o)^ I—*)**- 



V c -\^7t \ Jr h ••(«) V 



E 



|x3 = -/l 



\ 



(—he e ) u £ dx' 



J 
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Using (|34|) i . we have 

S ( D I 7fr^ dx ~ I 7h) Uedx 

V JY £ (uj) OX 3 OX 3 JY £ (uj) \ OX 3 J 

Dlim Y [ — (u-( x ',-h)-u+(x>0)) ut{x ^_ h)dx , 

e^o J r + {uA r° l / t(v eii )„\ 
ij e/« «r r °.«.« (w) / — jp- ; .--exp dt 

h {v e>ij ) a 



(u (x', -h) - u + (x', 0)) exp 

— Dlim > / — n -, ; (x 1 , 0) da?'. 

n^/W"l /" ? -L- a ^ exp f * 

Introducing the change of variables s = —£~ 6 t, we get 

1 *K,y) 3 



exp dt 



fc ft(-£-«t)j,(-£r«f) "V f 
= f k —L— exp f ds -+ /i (1/g 2 ) , 

using the ergodicity property (j2J) . Thus, using the proof of Lemma El we get 

lim | D [ — — ^ ^ u °' e dx — [ (fa, ^ U ° i u F dx 



V J Yc ( w )dx3 dx 3 " Jyjui) V £ ' S 3 dx 3 1 



t {u) ux 3 JY e {w) \ 

h (W) L r° {u ~ " M+) ~ Mo ~ " u+) exp ( 4V(vV) ^ 1 1 dx '- 

We now compute, using Lemma [13] 

lim I D I V T u £ ■ V T u , £ dx - / ((t> £ ) • V r w , e ) u £ dx J 

= lim ( / D*V T u £ ■ (V T u) 0e dx - / ( fa £ ) r • (V T w) 0£ ) M £ dx 

= h (q 2 (0)> / D*V T u+ ■ (V r «+) dx' - h (q 2 (0)> / ((v ) T • (V r u + )) <dx 
is is 



which leads to the limit variational formulation (1301) . ■ 

The problem associated to this limit variational formulation (130]) is given in the fol- 
lowing Corollary. 

Corollary 16 The sequence (u e ) E , where u e is the solution of r\-converges to the 
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solution uq of the problem 



- div (3 Vn+) + v+ d ■ Vn+ = \Z l \f mtt + , 
div (d\7uq ) + Vq A • Vmq =0 inVt h , 



-DVu + ■ e 3 
-/i (g 2 (0)) div r (D*Vi4) 

+ h (0)) K/)t • V,«f = ( ^- ( % „ ) on r } , (35, 

L>*V T n^ • n T = on 9E, 



+ 



h(l/q 2 (0)) V -Uu(.,-h)A 



on r+, 

= on r~, 



w 



w/jere A = exp ( g^f^y ^ j, «o |n+=: «o « nrf M o | n -= : M o ■ 

Proof. This is an immediate consequence of the limit variational formulation (I30I) . ■ 

Remark 17 Consider the case of a dispersive contaminant with a diffusion coefficient 
D (x, ui) defined through D (x, u) = D mo i + D disp (x, u) with 



Ddisp (x, u) 



D disp (x) in VL + U Q h , 

Ddisp (xi, x 2 , s~ d x 3 + atij (u) , u) in Y £)ij (u) , 



where (a^ - gZ is a sequence of random variables such that |ay (u>)\ < C, \/i,j G Z, 
w'^A probability 1, C being some non-random constant. We suppose that D is continuous 
with respect to the variable x and, with probability 1, do < D(x,u) < d\, where do and 
d\ are positive and non-random constants. We suppose that Ddi sp is a stationary random 
process. 

Let u G C 2 (O) be such that u = on T . We build the modified test-function u E) ij 
inside the fissure Y e ^ (co) 

U e ,ij (X!, X 2 , X 3 ) = U + (xi,X 2 , 0) 

(u~ (xi,x 2 , -h) - u + (x u x 2 ,0)) [ — _„ P \ — zfivit i ^fU\ dt 



+ 



•T3 



,qi{-e-H) qj {-e-H) D^-e'H) 



exp (-t (v £ ,ij) 3 ) 

_ hqi (-e-H)q 3 {-e-H)D*{-e-H) 

where (D*) (s) = {D mo i + D disp (ie,je, s + (u) , u)). Implementing this test-function in 
the above process, one gets at the limit a problem similar to (T23J), except that (1/q 2 (0)) 
is now replaced by (.) q 2 (0)), where (1/ (D^q 2 (0))) is the mathematical expectation 
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of 1/ (q 2 (t) (-,t)), with respect to the measure probability P, and bj is replaced by the 
Z -periodic solution bj of the problem 

div (D (x) (e-j + Vbj)) = in Z 1 j = 1, 2, 3, 
(V6j + ej) ■ n = on S, 

and c m by 

div (D (x, lo) (e m + Vc m )) = in Z' m = 1,2, 
(Vc m + e m ) ■ n = on dZ'. 

4.3 The asymptotic behaviour in the case of a reactive contam- 
inant (11 > 0) 

In this subsection, we consider the reaction-diffusion equation §5§ with first-order reaction, 
that is with 1Z > (see for example [3], [II]). We denote w £ti j the solution of the 
differential equation 

~ D i~ 3 ( (<Mj) (~ rts ) ilf ) (X3) ~ ((Mj ) (- e ~ fla *) ^ilf (xs) 

+TZ (qiqj) (-e~ e x 3 ) w £:ij (x 3 ) = 0, 

<o-(o) = o 

(36) 

and the solution of 

<9 / <9z ' ■ \ dz ■ ■ 

+^ {QiQj) (-£~ x 3 ) z £!ij (x 3 ) = 0, 
* e ,«(0) = 

<«(0) : 



1 



where is the velocity defined in (l3~Tj) . We have the following estimates. 



QiQj (°) ' 
(37) 



Proposition 18 There exist non-random positive constants Co and C\ independent of e 
and of i and j , such that: 

1. Me > 0, V(z,j) e I e (u), Vx 3 e [-h,0] : 1 < w e>ij (x 3 ) < C x and -C^x 3 < 
Z £ ,ij (x 3 ) < d. 
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2. V{i,j)el e (u),Vx 3 e[-h,0] 



lim 

£^0 



lim 

e-i-0 



lim 



K,ij ( x s) (<Mj) (-£ e x 3 ) exp ( 



z s,ij \%3 
•^3^e,ij 



lim 

£->0 



V D 

z 'e,ij M (fcty) (-£~ 6 x 3 ) exp 



w e ,y (x 3 ) - cosh [Tlxs 

n v 

sinh ( lZx 3 



D 



— cosh TZx-s 



0. 
0. 



where K = ^TZ{q 2 (0)) (1/q 2 (0)) /D. 



Proof. 1. Multiplying the equations ( 155|) and ( 157)) by exp (x 3 v Eii j/D) and integrating by 
parts the first term of these equations, we obtain the Voltera type integral equations 



W e,ij [X3 ) 




(QiQj) (~ £ 6s ) ex P ( s ^d~ 



exp -C 



D 



d( \ds + 1, 



(QiQj) (~ £ ~ e 
(J (QiQj) {~£~ 9 s) ex P (s^ 




exp -C 




(3f 



23 exp 



(QiQj) (-e- e C) 



(which can be solved by the method of successive approximations). Taking into account 
the hypothesis (CQ), we obtain the first point of the Proposition. 
2. Consider the integral equations 



K 



v>ijM = ^(? 2 (0))(l/g 2 (0)) / (x 3 -s)w ij (s)ds + l 



D 
K 



*vM = ^(g 2 (0))(l/g 2 (0)) / (x 3 -s)z tJ (s)ds + x 3 (l/q 2 (0)) 



whose solutions are Wij (x 3 ) = cosh yZx 3 j and Zij (x 3 ) = (1/q 2 (0)) sinh \lZx 3 J /TZ re- 
spectively. The differences w e ^j — Wij and z £ ^j — satisfy integral equations of Voltera 
type. Using the ergodic result ©, we prove the convergences 

lim \w e>ij (x 3 ) - w i:j (x 3 )\ = lim \z £>ij (x 3 ) - z {j (x 3 )\ = 0, 



e^0 



e-s-0 



uniformly with respect to x 3 E [—h, 0]. 

The last estimates for the derivatives follow from the derivation of the equations (138)) . 
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Let u G C 2 (JTj be such that u = on T + U T . We here define the test-function 



inside the fissure Y E) ij (lo) through 

U £ ,ij (x) = U + (Xi, X 2 , 0) W £tij (x 3 ) + 



u (xi,x 2 , —h) - u + (xi, x 2 , 0) w e4j (-h) 

Z e ,ij (~h) 



We finally define the test-function u 0j£ in the same way as (133]) . It is easily proved 
that the sequence (uq )E ) "/^-converges to u. On the other hand, we compute 

u 0i£ (x',0) = u + (x 1 ,x 2 ,0), 
u £ ,ij (x', -h) = vr (x u x 2 , -h) , 
9«o, £ , u~(x u x 2 ,-h) -u + (x 1 ,x 2 ,0)w e>ij (-h) 1 



o 2:3=0 

dx 3 



du . 



z e,ij (~ h ) 



QiQj (0) ' 



-q^- \ X3 =-h = u + (x 1 ,x 2 ,0)w' £ ij (-h) 



| u (x u x 2 , -h) - u + (x 1 ,x 2 ,0)w eiij (-h) , , 



z £ ,ij (~h) 



We have the following result. 
Lemma 19 One has: 

dUn £ 



lim > / 

tl< Jri 



,. , , Jr+ ..(u) 9x 3 3 



=0 (ffity) (0)tt+do- 



-1 



lim > / 



V? 2 (0))7s sin h (ftA 

<9w , £ I 



w — ?i + cosh ( IZh ) ) Uq(1x', 



V -\^7t \ J r~ ,.(u) 9x 3 3 



1 



K 



(1/q 2 (0))7 2sinh 



-ft (QiQj) (he 9 ) u £ da 

u~ cosh (iZh) — n + ) exp (777-7^77 — ) u^dx', 



D(l/q 2 } fi 



Proof. Using the estimates of Proposition [181 we can replace, when e is small enough, 
w e,ij (—h) by cosh [TZh), z £i j (—h) by — (1/q 2 (0)) sinh (7£/i) /7£, w' eij (—h) by the quan- 



tity 



exp(hv e ,jj/D) Tl 
' {q iqj )(e-<>h) <l/9 2 (0)> 



sinh [TZh] and 2^. (—h) by the quantity 



exp(hv e! jj/D) 
i<li<lo){ e ~ eh ) 



cosh f ) . 
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Thus 



du 0j£ 
dx 3 



\x 3 =-h (qiQj) {he e ) 

u + (x ll x 2 ,0)w' £jij (-h) \ 
i vT (x 1: x 2 , -h) - u + (xi, x 2 , 0) w £ , tj (-h) , 



£^0 



\ 



-K- 



{m){e-°h) (l/g 2 (0)) 
it~ (xi, x 2 , —h) — m + (xi, x 2 , 0) cosh ^7£/ij 



X 



{1/q 2 (0)) sinh ( Kh 
exp (hv £>ij /D) 



\ ' {Wj){e- 6 h) 



K 



exp (hv £> ij/D) 



(m) (he- 9 ) 



^ — u + (xi, x 2 , 0) sinh (jZh^j \ 
u~ (xi, x 2 , —h) — u + (xi, x 2 , 0) cosh ( IZh 



cosh IZh 



sinh f IZh 



K 



y x cosh (iZh^j 
exp (hv £:ij /D) 



(1/q 2 (0)) sinh (jih) 

l + (xi,x 2 ,0) (- (sinh) 2 (ft/i) + (cosh) 2 (^)) 



— -u (xi, x 2 , — h) cosh ( 

n 



(1/q 2 (0))sinh (ft/i 



exp (hv e jj/D) (u + (xi,x 2 , 0) — u (x\ ) x 2) —h) cosh ^7£/i^ , 



whence 



du 0>e 



lim > 

e->0 



I .'C'l 



-1 



<W(o)>y Ssinh (nh) 



=-h (qiQj) {he e ) u £ da 
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Similarly, we get 



Jr+ (<„\ 9x3 



lim > 



I a^3 : 



=0 {m 3 ) (0) u+da 



lim > / 



(ij')e/ s (w) ^o^.ii 



lim > 



u~ (xi, x 2 , — /i) - n + (xi, x 2 , 0) w £>i j (—h) 1 
U~ (.Ti, x 2 , -/i) - u + (Xi, x 2 , 0) w sij (-ti) 



(QiQj) (0) n+da 



Jr+ a 



z e,ij y-h) 

u~ (xx, x 2 , —h) — u + (xi, x 2 , 0) cosh ( TZh 



-1 



TZ 



(W (0))7s sin h (Kh 



{1/q 2 (0)) sinh [TZh 
u~ — u + cosh [TZh] ) utdx' 



■ufda 



which leads to the desired limit. ■ 

Now, using the same methods as in the above subsection, we obtain the following 
Theorem. 

Theorem 20 The sequence {u £ ) e , where u e is the solution of |3]], Tx-converges to the 
solution u of the problem 



- div (3Vu+) + v+ d ■ V< + TZu+ 

- div ( D\7u ) + v o d ■ Vn + TZuq 
-DVu+.e 3 - h (q 2 (0)> div T (D*Vu+) 

+h(q 2 (0)> («b f/ ) T • V T «+ 



\Z l \f 




( D 
\ 



in Q + 
in Q 



h > 



TZ 



\ 



W (0)) sinh (Oh 



x ( u o cosh ( TZh ) — u A I y 



on r 



o J 



DVu ■ e 3 







( ° 
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<W(0)> sinh TZh 



x ( n cosh ( TZh ) A — Mit 



on <9E 
on r 








on r + , 
on r~, 

(39) 



wifli A = exp ( gffi)^) ^), Uq ■= uo |n+ and n := n | n - 
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